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Modeling of a heat equation with a Dirac density 


Scott W Hansen * and Jose de Jesus Martinez * 


Abstract 

We consider a linear hybrid system consisting of two rods connected by a thin wall of 
width 2e and density l/2e. By passing to a limit, we obtain a system describing heat flow 
of two rods connected by a “point mass” whose dynamics are governed by a differential 
equation. We prove that solutions of the epsilon problem converge weakly to solutions of 
the limiting system. 


1 INTRODUCTION 


Consider a linear hybrid system consisting of two wires or rods connected by a thin wall of 
width 2e > 0 and density l/2e. Assume the two rods occupy the intervals = (—Ai, —e) and 
a;e _2 = (e, L 2 ), and the wall occupies the interval We = (—e, e). Correspondingly, let x), 

Ve = v^(t,x) and = Ze{t,x) denote the temperature distribution on their respective domains 
and Wf. We suppose the temperature of the rods and wall satisfy the heat equation 
on their respective domains with Dirichlet boundary conditions at endpoints x = —Li^ L 2 . The 
linear equation modeling heat flow of such a system is as follows: 


Cipilie — kiUe" = 0 , 

C2P2Ve - k2Ve" = 0 , 

- kze" = 0 , 

u^{t, -e) = Zeit, -e), Ze{t, e) = v^{t, e), 
kiu^'{t, -e) = kze{t, -e), e) = e), 

Ue{t, -Li) = Ve{t,L2) = 0, 


t > 0, X S Wep 
t > 0, X G Wep 
t > 0, X G UJe 
t > 0 
t > 0 
t > 0. 


( 1 ) 


Throughout this article, ' will denote spatial derivatives and ' will denote temporal derivatives. 
The parameters c > 0 and A: > 0 in the third equation represent the specific heat and con¬ 
ductivity of the wall connecting the two rods. The parameters cp pi and ki in © are positive 
and represent the specific heat, density and thermal conductivity of the rod on the subdomain 
Wep. It will later be convenient to use the diffusivity coefficient af = kiJciPi for i = 1,2. The 
fourth equation guarantees continuity of the temperature across the interface x = ±e and the 
fifth equation represents the heat flux continuity condition at the interfaces (see m Chapter 
8] ). We complete the system by adding the initial conditions 

{u°(x),u°(x),z°(x)} = {uq0,x),'uq0,x),zq0,x)} (2) 


in an appropriately defined function space at time t = 0 so we may determine the solution of 
o uniquely. 
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We show in this article that with appropriate assumptions on the initial conditions, the 
solution {Ue,Ve,Ze} of (HD with ([2]) converges in a weak sense to the solution of the following 
limiting hybrid system; 


cipiii — kiu" = 0, 

C 2 P 2 V — k2v" = 0, 

< c i = k 2 v'{t, 0 ) — 0), 

u(t, 0) = v{t, 0) = z{t), 
u{t, -Li) = v(t,L 2 ) = 0, 


t > 0, 
t > 0, 
t > 0 
t > 0 
t > 0, 


X G OJi := Wo,i 
X G 0J2 ■= Wo,2 


with initial conditions of the form 


( 3 ) 


{u°{x),v°{x),z°} = {M(0,a;),u(0,a;),z(0)} 


( 4 ) 


given in an appropriately defined function space at time t = 0. The third equation in ([31 states 
that the rate of change in temperature of the point mass is proportional to the net heat flux 
into the point mass. This can be viewed as a form of Tick’s law of diffusion. 

Similar hybrid systems involving strings and beams with point masses have been studied 
in the context of controllability and stabilization theory. See for example 0, m, m, m, uni, 
and 0. In particular, C. Castro showed in 0 that a system similar to m 
with strings can be obtained from a system similar to that in (HD and gave a detailed spectral 
analysis. 


2 WELL-POSEDNESS 

We begin by proving well-posedness of the limit problem ([3]) . 


2.1 The limit problem 

Given u, v defined on wi, W 2 and 2 ; G R, let y = (u,v,zy, where * denotes transposition and 
define 


n = L^iuji) X L^{u}2) X K 


equipped with the norm 


\\y\\n = ClPlMl^ + C2P2\\v\\l^ + c Izp 
where || • is the usual norm in LF‘{u!i) for i = 1,2. Define 

■dull = {m e I Zi(-Li) = 0} 

■0^^ = {U G H^{ijJ2) I v{L2) = 0} 
d = {(u,!)) G d^^ X I u(0) = z)(0)} 


equipped with the norms 



ki\\u' 


2 

L^(uii)i 


l^ = N 


Iklll 


for z = 1,2. We can check that (see 0) is algebraically and topologically equivalent to Hq{Q), 
however one can think of as a subspace of d^n x daj 2 ■ The space 


W = {(it, z;, z) G X R I u(0) = u(0) = z} 
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is a closed subspace of i? x R with norm we may define as \\y\\w = IK'^jI t is easy to see 
that the space W is densely and continuously embedded in the space %. Define the unbounded 
operator A : D{A) C "H —>■ "H by 


/ al(P 0 0 \ 

= j 0 a2(f 0 I (5) 

\-^5od ^Sod 0 / 

where d denotes the (distributional) derivative operator and Sq denotes the Dirac delta function 
with mass at a; = 0, and the domain D{A) of A is given by 

D{A) = {y€W : u€ i?^(wi), v € i?^(w 2 )}. (6) 

When D{A) is endowed with the graph-norm topology 

\\y\\liA) = Mn + Uyrn 

it becomes a Hilbert space with continuous embedding in 'H. We can therefore write the limit 
system m as 


y{t) = Ay{t), y(0) =y°, t> 0 (7) 

where y° = Since D{A) is dense in W and the latter is dense in "H, it follows that 

is a densely defined operator. 

Lemma 2.1.1. The operator A : D{A) TL is a bijection. 

Proof. Let / = (/, g,h) € H be arbitrary. Then the solution to Ay = / is given by 

/Cuix + Li) - F{x)\ 

y=\c,(x-L 2 )-G{x)\ ( 8 ) 

where 

-chL2 + k2{F{Q)-G{Q)) 
k2Li -\- kiL2 
chLi + fci(f(0)-G(0)) 

^2^1 -|- kiL2 

chLiL2 L2kiF{{f) L\k2G(Qi) 

^2^1 -t- k\L2 

Since {u",v") = S L^(wi) x L^(w 2 ) it follows that {u,v) S x H'^{uj 2 )- 

Furthermore, one can check from that u(0) = u(0) = z and u{—Li) = v{L 2 ) = 0 so that 
y G D{A). Thus A : D{A) —?> H is surjective. 

Finally, note that the null space of A is trivial since when / = (0, 0,0) we see that y is the 
trivial solution. Then A is injective and hence bijective. □ 

Lemma 2.1.2. The operator A : D[A) % is symmetric and dissipative. 



px pO 

Gu 

F{x) = 

/ / ar^/(^) drds, 



J-LiJs 

pL2 ps 

Gx 

G{x) = 

/ / drds, 



Jx Jo 

a 
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Proof. Consider (p = (^, z/, G D{A). Then 

{Ay,(p)-H = A:imVI°Li - ki{u',p')u,i + k 2 v'iy\Q^ - k 2 {v',A)^^ + {k2v'{0) - kiu'{0))C 


= -kiup'\_i^^ + ki{u,pL")ui^ - k2vv'\^^ + k2{v,v") 


= cipi(u,af/r")^j + C2P2(^',Q:2^^")^^2 + cz ( - —m'(0) 

V c c 

= {y:-^p)n 

for all y = {u,v,z) G D{A). Hence D{A) C D{A*) and so ^ is a symmetric operator. In 
particular, when we choose y = p we see from the above computation that 

{Ay,y)n = -k4ur.,-k2\W\\l, = -\\yry,. 


^2 


Thus we have that {Ay, y)u < 0 for any y G D{A) as needed to show A is dissipative. 


□ 


Lemma 2.1.3. The operator A is closed, self-adjoint and its inverse is a compact operator in 

n. 

Proof. As mentioned before, A is densely defined in TL and from Lemmas 12 . 1.21 and 12 . 1.11 we 
have that it is symmetric and R{A) = TL. It follows from Theorem 13.11 in [T5], that A is 
self-adjoint and its inverse A~^ is bounded in TL. Furthermore, since the inverse is bounded, 
we have 0 G p{A) and Theorem 13.9 in [15] implies that A is closed. 

Next we claim that K := A~^ is compact. From formulas (I8|)-(l9|) we can decompose 
K = Ki -f K 2 where 

f-F{x)\ /Cu{x + Li)\ 

Kif = I -G(a:) 1 , K2f= \ L 2 ) j ■ 

Since the mappings / !->■ F{x) and g !->■ G{x) are Volterra-type operators, Ki is compact. Since 
0 G p{A) and Ki is compact, K 2 must be bounded. Since also K 2 has finite rank, it follows 
that K 2 is compact, and hence also K is compact. □ 


Proposition 2.1.1. The operator A is the infinitesimal generator of a strongly continuous 
Co-semigroup of contractions which extends for Re{t) > 0 to an analytic semigroup. 

Proof. By Lemma [2.1.31 we have A is a closed, densely defined and self-adjoint operator and 
by Lemma 12.1.21 we see that both A and A* are dissipative. Therefore, by the Liimer-Phillips 
theorem (see Luo et al [Sj), we have that A generates a Co-semigroup of contractions. 

Furthermore, the spectrum o'(A) of A is contained in (—c»,0) and by the computations 
shown in m page 55], we obtain ||i?(z,A)|| < sec(0/2)/|z| for all z = pe*® G C \ (—oo,0] and 
9 G (7r/2, tt). Rewriting the angle 9 as tt/2 -\- 6 where 0 < <5 < tt/2 and letting M = sec{9/2) we 
see that 

11 ^( 2 :, A)II < ^ for all z €Sg 

FI 

where 


Sg = {z G C : |argz|<d}. 

Note that Sg U {0} is contained in the resolvent set p{A) of A. By Theorem 5.2. in [Tl] we 
have that T can be extended to an analytic semigroup in a sector 85 . If Re(z) > 0 then for 
6 < t:/ 2 large enough, that z G 85 . Hence A generates an analytic semigroup in the right half 
plane Re(z) >0. □ 


4 










As a consequence of Proposition 12.1.11 given initial data G there exists a unique 
solution 


2 /eC'([ 0 ,oo);?^) (10) 

to the Cauchy problem ([T]). If in addition, ?/° G D{A), then y G C([0,oo)- D{A)). 

2.2 The approximate problem 

Now consider functions u^, Ve and defined on uje,i, uje ,2 and Wg respectively and define 
= {ue,Ve,Zey. Let 

Tie = x L\u;,,2) x L^uj,). 

equipped with the norm 

+ C2P2\\vy\l^, + 

where || • . is the usual norm in L?{ujc,i) for z = 1, 2 and || • is the usual norm in L^{uje). 

Define 


= {Ue G I Ue{-Li) = 0} 

e H^ii^e,2) I VeiL2) = 0} 


equipped with the norms 

for i = 1,2. Next, consider the following subspace of ^ x ^ x 

We = {ye G X 2 X H^{uJe) I Uei-e) = Ze{-e), Ze{e) = Ue(e)} 


with the norm 


l2/.llk=fcillV|ll,,+fc2||V|ll., 


/||2 


Remark 2.2.1. It is easy to show that the spaces We are uniformly equivalent to in 

the sense that there exists some constant C > 0 such that 


^llv^llw. < ll‘dllffi(n) < C'll'dilw, 

where C is independent of e for all 0 < e < eq with finite eg. Furthermore, it is easy to see that 
the space We is densely and continuously embedded in the space "He. 

We will also make use of the space = H‘^{u}e,i) xH'^iuJe, 2 ) xH‘^{ijJe)- Define the unbounded 
operators Ae ■ D{Ae) C "He —?> 'He by 

[alcP 
Ae=\ 0 

V 0 


0 0 \ 

al<P 0 , ( 11 ) 

0 
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with domain D[Ae) given by 

£)(A) = [ve G We : ye G "He: kiuj{-t) = fcZe'(-e), fcZe'(e) = fc2t’e'(e)} ■ 

When D{Ae) is equipped with the graph-norm topology 

lly.llD(^.) = lly.llH. + IIAye||?,,, 

it becomes a Hilbert space with continuous embedding in "He- We can now rewrite system (ED 
as a Cauchy problem: 

ye{t) = Aeye{t), ye(0) = G "He, t > 0. (12) 

It is easy to see that A^ is densely defined on As in Section I^TT] we have the following 
results. 

Lemma 2.2.1. The operator Ae ■ D{Ae) Tie is a bijective, dissipative, closed, self-adjoint 
operator with a compact inverse in He • 

Proposition 2.2.1. The operator Ae is the infinitesimal generator of a strongly continuous 
Co-semigroup of contractions which extends for Re(t) > Q to an analytic semigroup. 

The fact that Ae is the infinitesimal generator of an analytic Co-semigroup implies that for 
all y'^ G He there exists a unique solution 

y, GC([0,oo);iH,) (13) 

to m- Moreover, if y° S D{Ae), then also ye G C([0, oo); Il(v4e))- 


3 WEAK CONVERGENCE 


The energy functional of the hybrid system m is given by E{t) = ||y||^/2. By taking test 
functions ip G Co([0,oo) x fl), a weak form of the hybrid system (|3]) is given by 

J Cipiu^ip{0,x) dx-\- J C 2 P 2 V^p( 0 , x) dx -\- cz^(p{0,0) 

UJl i02 

OO 

= —y|y cipiwp dx + J C 2 P 2 V(p dx-\-cz(p{t,0)'^ dt ^^4) 


0 u;i 


kiu'(p' dx-\- k 2 v'p' dx > dt. 


On the other hand, the energy functional for the e-dependent problem ED is Ee{t) = \\yernj2 
and by taking test functions p G Cq ([0, oo) x O), a weak form is 

y cipiu^p{0,x) dx-\- y C 2 P 2 VePi 0 ,x) dx-\- y -^z°p{0,x) dx 




0 

OO 


-b 


0 LJ,,! 




CiPiUeP dx -\- J C2P2VeP dx + J dx ^ dt 

klUe'p' dx-\- J k2Ve'p' dx-\- J kze'p' dx )■ dt. 

<^e,2 <^€ 


(15) 
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We give sufficient conditions such that we may pass to the limit in (TTSl) to consequently obtain 
m- Assume that G D{Ae) and furthermore, there exists Mi > 0 such that 

\\y°\\n.<Mi, (16) 


for all e > 0. Then we obtain the following result. 

Lemma 3.0.2. The energy of system 0 is (uniformly) bounded by the initial energy in the 
sense that there exists constant C > 0 such that E^ft) < C for all e > 0 whenever the initial 
data y) satisfies m- 

Proof. Note that if y) G D{Ae), then G C'([0, oo); D{At)) and the energy satisfies 


Mt) 


+ k2Vc'Ve\^^ + 

— j ki\uf'\^ dx — J fc 2 |Ue'P dx — J 

l^e.l <^€,2 

-lly.llk. 


fcjZe'P dx 


(17) 


which implies i?e(t) < 0, and thus Ef^{t) < Ec(0). Hence by density, there exists C = Mf/2 
for which £'e(t) < C for alH > 0 and initial data satisfying (ITSl) . Consequently, we find that 
ye G L°°([0,oo);'H,) for all e > 0. □ 


Now assume there exists M 2 > 0 such that 

IlyellWe < M2, ( 18 ) 

for all e > 0. Then we obtain the following result. 

Lemma 3.0.3. Assuming condition (1181) holds, the sequence solutions {ye}e>o to problem (HD 
is uniformly bounded in L°° ([ 0 , 00 ); 77^(11)). 

Proof. Since G D{Ae) we have that G (^([0, 00 ); Il(Ae)) and 
^^lll/ellw, = k2V^vf\^^'^ + kZeZe'[_^ 

- fci(?ie,Ue")e,i dx - fc2 (u'eVe")e,2 “ k{ZeZe'')e 

= -|| 2 //|lk- 

This shows that the sequence {||ye(f)||w }t>o is monotone decreasing in t and thus by density, 
||y£(f)||we < ll?/ellwe < M 2 for all t > 0 as needed to show y^ G L°“(0, 00 ; We). From Re¬ 
mark [2T2T] we see that the spaces We and i7o(^) equivalent and thus there exists K > 0 
independent of e such that 

ll2/e|lL“([0,oo);ffRn)) ^ ^ 

for all e > 0, as needed to show solutions to the e-dependent problem are uniformly bounded 
in 3 ;. □ 


Next, it is natural to assume the initial data convergences in a weak sense in Tie. It is 
easy to see that being a sequence in L^{uje,i) implies that is a sequence in Lffwi). 

Likewise, Xuj^. 2 ' 0 ^ ^ L^(wi). Then we will assume that 

{ Xwe.Wc ^ '“0 weakly in Lf{uJi) as e —>■ 0 

Xwe, 2 ^e ^ ■^o weakly in L‘^{uj 2 ) as e —>■ 0 (19) 

^ dx —>■ zq in IR as e —>■ 0 . 

The following theorem is our main result. 
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Theorem 3.0.1. Let {2/e}£>o be the sequence of solutions to the e—dependent problem o with 
initial data y^. Assuming (HU), dni) and (ITOl) . the family {yf}^yo of solutions to (HD problem 
satisfies 


ye ^ y in L°°{[0,oo-,Hl{n)) 

as e 0 where y is the weak solution to the limit problem ([3]) with initial y’^ G TL. 

Proof. From Lemmas 13.0.21 and 13.0.31 the initial energy provides a uniform bound for the 
solutions ye G L°°([0, oo; Lfg (^)) ■ We can then extract a subsequence of solutions (which is 
still denoted by the index e) such that 

u in L°^(0,oo;L^(cui)) nL°°(0,oo;-de^J (20) 

^ u in L°“( 0 ,oo;L^(a; 2 )) nL°°(0, oo;d,.„ 2 ). (21) 

Next, observe that ge{t) := ^{l,Ze{t))e defines a function on [0,oo). Applying Holder’s in¬ 
equality we see that \ge{t)\ < ||ze||e/-\/^- By condition (TTHl) we have that {g£}e>o is a uniformly 
bounded sequence in L°°{0, oo). Invoking the Banach-Alaoglu Theorem we can extract a sub¬ 
sequence of Z£ (still denoted with the index e) and find z G L°°{0,oo) such that 

ge zinL°°{0,oo) ( 22 ) 

as e —>■ 0. We now pass to the limit in each of the nine terms in the characterization (fT5)l of 
weak solutions of the e-problem with ip G (^(^([O, oo) x H). Since (/3(0, •) G Cq)!!) it follows from 
assumption HID on the initial data y^ that 


/ u^ip(0,x) dx= dx^ u^ip{0,x) dx, 

OJe,! J J UJl 

/ u°(/3(0,a;) dx= 2U°(/3(0, cc) dx ^ i v^(p(0,x) dx. 

JcJe ,2 ^ OJ2 ^ OJ2 

Next we claim that ^z^(p(0,x) dx —>■ z°(/5(0,0). By adding and subtracting the term 
z^ip{Q, x) and applying the triangle inequality we find that 


/ “ ^V(0, 0) 


< max|(/?(0,a;)| 


X^UJ, 


Z° - z° dx + \z°\y J |</5(0,a;) - (/7(0,0)| dx. 


The first term in the right hand side above tends to zero from HID and the boundedness of 
(/7(0, •) on uje- The second term tends to zero as well by the continuity of (/^(O, x). 

From (1301) and (1311) we have that in particular for (p G C{\f), 00 ) x H) that 


0 J iiJ 


UeP dx dt - 
VeP dx dt 


/ 0 J uji 


u(p{t, x) dx dt, 


JQ 

Next we want to show that 


/o ^ UJ2 


v(p{t, x) dx dt. 


1*00 ^ I * 1*00 

/ — Ze{t,x)(f{t,x) dx dt ^ / z{t)(p{t,0) dt. 

Jo 2e Jq 


(23) 
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Observe that 


/o 


oo 


POO 

z^{t,x)(p{t,x) dx dt — / z{t)ip{t,0) dx dt 

Jo 

f°° 1 f 

/ — Zeip{t,x) — Z(p{t,0) dx dt 

Jo 2e 

Ze{(p{t, x) — ip{t, 0)) dx dt 


1 




(zg — z)(fi{t, 0) dx dt 


The last term tends to zero by (|22]) . Regarding the first term, observe that applying Holder’s 
inequality we have 

[ ^ [ z^(ip{t,x) - (p(t,0)) dx dt < f ^ f |zg| - (/)(t,0)| dx dt 

Jo Jui^ Jo ze 

~ Iq - (^(t,0)| dx dt 


< 


\/2i 


11^.11. 




\ip(t,x) — (p{t,0)\ dx dt. 


Note that || 2 ;g||e/-\/ 2 e is bounded by condition (IT6t and by the continuity of (p, we have that the 
above tends to zero as e —>■ 0. Since p' € C([0,oo) x H) we have from (l20)l and (I2T11 that 


^OO 


^OO 


/ / Ue (p dxdt= / 

'0 Jq Juji 

poo p pOO p 


Xiu, ^ujp' dx dt 


u'p' dx dt 


'0 ^ LUl 


vjp' dxdt= / Xuie dt 


v'p' dx dt. 


' 0 J u), 


€,2 


/O J (jJ2 


/o j U )2 


Finally, the term kzjp' dxdt tends to zero as e —?> 0 since Ze{t,-) g and p' £ 

C([0,oo) X H). 

From the above discussion we now have that there exists a convergent subsequence of 
in the weak star sense, whose limit satisfies the equation da. Since the limiting system has 
a unique weak solution, it follows that the convergence holds for the whole sequence {ye}e>o- 
Therefore, we have shown that the limiting system © can be approximated with the sequence 
of e dependent problems H]) as needed. □ 
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